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Complexity in materials often arises from competing interactions at the atomic length 
scale. One such example are the strongly correlated heavy-fermion materials where 
the competition^ between Kondo screening^ and antiferromagnetic ordering is believed 
to be the origin of their puzzling non-Fermi-liquid properties^'- . Insight into such 
complex physical behavior in strongly correlated electron systems can be gained by 
impurity 

dopin^HUl. 

Here, we develop a microscopic theoretical framework to demon- 
strate that defects implanted in heavy-fermion materials provide an opportunity for 
unveiling competing interactions and their correlations in real space. Defect-induced 
perturbations in the electronic and magnetic correlations possess characteristically 
different spatial patterns that can be visualized via their spectroscopic signatures in 
the local density of states or non-local spin susceptibility. These real space patterns 
provide insight into the complex electronic structure of heavy-fermion materials, the 
light or heavy character of the perturbed states, and the hybridization between them. 
The strongly correlated nature of these materials also manifests itself in highly non- 
linear quantum interference effects between defects that can drive the system through 
a first-order phase transition to a novel inhomogeneous ground state. 



The essential features of heavy fermion materials is 
a Kondo lattice of magnetic atoms with localized mo- 
ments arising from 4/ or 5/-electron shells, and a delo- 
calized conduction bancP^. The magnetic moments are 
coupled to each other via an antiferromagnetic interac- 
tion, J > 0, and to the spin moments o f th e conduction 
electrons via the Kondo coupling, J > Cf^^. The ground 
state of this system is realized either by Kondo screening 
of the magnetic moments, resulting in a heavy-Fermi- 
liquid phas^^, or by their antiferromagnetic ordering. 
These two competing phases, which manifest strongly 
correlated many-body states, are separated by a quantum 
critical point ( QCP). To date, no theoretical consensus 
has been reachecP^H^on the physical origin for the aston- 
ishing deviations from Landau's Fermi-liquid theory ob- 
served in the quantum critical region above the 
Understanding how the formation of antiferromagnetism 
competes locally with the creation of a Kondo singlet, 
how coherence is established in Kondo lattices, and how 
the resulting magnetic and electronic correlations are in- 
tertwined, is crucial in unraveling this mystery. 

By their very nature defects implanted in a Kondo 
lattic^^nmH jjianifest local perturbations that can reveal 
the form of strong correlations in real space. These de- 
fects can take the form of missing magnetic Kondo atoms, 
i.e., Kondo holes, or of non-magnetic impurities. By dif- 
ferently perturbing the electronic and magnetic states of 
the system, defects provide an unprecedented opportu- 
nity to differentiate (in real space) between electronic cor- 
relations arising from Kondo screening, and the antifer- 
romagnetic correlations between the magnetic moments, 
as schematically shown in Fig. SI of the Supplementary 
Information. The spatial form and extent of these per- 
turbations are not only a measure for the strength of the 
correlations, but also provide insight into the hybridiza- 
tion of the (heavy) magnetic /-electron states with the 



(light) conduction band. This hybridization results in 
the Kondo screening of the magnetic moments and the 
opening of a gap in the electronic spectrum. The experi- 
mentally accessible, spectroscopic fingerprints of the per- 
turbed electronic and magnetic correlations can be found 
in a redistribution of spectral weight in the local density 
of states of the conduction band, and in an enhancement 
or suppression of the /-electron spin susceptibility, re- 
spectively. Modulations in the density of the conduction 
electrons, as induced by a non-magnetic impurity, have 
profound effects on the spatial form of the correlations, 
and lead to the formation of an impurity bound inside the 
hybridization gap. The coupling between electronic and 
magnetic correlations can create highly non-linear feed- 
back effects on quantum interference in periodic arrays 
of Kondo holes, leading to the spatial reconstruction of 
the ground state via a first order phase transition. Our 
findings open a new path for studying competing interac- 
tions in real space, and provide a microscopic, real space 
picture of the perturbations induced by defects in heavy 
fermion materials. 

In order to study the spatial perturbations induced by 
defects in the magnetic and electronic correlations of a 
(heavy fermion) Kondo lattice system, we have developed 
a real-space formalism, the Kondo-Bogoliubov-de Gennes 
method, whose details are described in Sec. S2 of the Sup- 
plementary Information. This theoretical formalism pro- 
vides us with two salient quantities that determine the 
physical properties of the system. The first one is the 
site-dependent hybridization, s(r), between the conduc- 
tion band and the /-electron state of a Kondo atom. It 
is a measure of the Kondo screening of the magnetic mo- 
ment, and of the electronic correlations in the resulting 
many-body state. Here, we focus on the Kondo screened 
phase of the heavy-fermion materials where s(r) ^0 for 
all sites. The second quantity, x(r,r')? represents the 
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Figure 1: A hole in the Kondo lattice a. Schematic 
picture of a Kondo lattice containing a Kondo hole: green 
(orange) dots represent the sites of the conduction lattice 
(magnetic atoms), b. Electronic Structure of an unperturbed 
Kondo lattice system for I / J = 0.001: the large Fermi sur- 
face (black) arises from the hybridization of the magnetic /- 
electron states with the conduction band, the small Fermi 
surface (red) represents the (unhybridized) conduction band. 
Spatial contour plots of c. As and d. Aric for I / J — 0.001. 
Spatial contour plots of Ax for e. I/J = 0.0005, and f. 
I/J — 0.0011. The relative changes, As, Anc, and Ax, are 
obtained by subtracting the results for s(r), nc(r) and x(r, r') 
in the unperturbed Kondo lattice from those of the Kondo 
hole system, and dividing by the former. We display the 
value of x(r, r') between nearest neighbor sites r and r' at 
(r + r')/2. 



strength of short-range antiferromagnetic correlations be- 
tween inagnetie moments on (nearest-neighbor) sites r 
and r<™ 

The simplest defect in a Kondo lattice is realized by 
removing a magnetic Kondo atom from it, resulting in 
a Kondo hole, as shown in Fig. [T^. To demonstrate the 
spatial perturbations in the electronic and magnetic cor- 
relations of the system that are induced by the Kondo 
hole, we present in Figs, [ij; and [T}i, a two-dimensional 
contour plot of the relative change in the hybridization, 
As, and the conduction electron density, Anc, respec- 
tively, between the Kondo lattice with and without a 
hole. Both quantities exhibit similar spatial oscillations, 
whose isotropy and wavelength of Sag (ao being the lat- 
tice constant) imply that they are determined by the 



Fermi surface of the unhybridized conduction band (see 
Fig. [T|3) with Fermi wavelength = lOao. These oscil- 
lations decay exponentially but remain practically un- 
changed over the range of I/J considered here. The 
spatial plots of Ax in Figs, [ijj and [if reveal strongly 
anisotropic oscillations along the lattice diagonal. Their 
weaker reflection can also be found in As (see Fig. [l|;), 
clearly demonstrating the coupling between the system's 
electronic and magnetic correlations. The spatial pattern 
of Ax is determined by the strongly anisotropic Fermi 
surface of the hybridized system, shown in Fig.[T]3, which 
possesses a large degree of nesting and a Fermi velocity, 
vi, along the lattice diagonal which is about 10 times 
larger than that along the bond direction, V2- This con- 
clusion receives further support from the rapid oscilla- 
tions in Ax along the lattice diagonal (see Fig. S2 in the 
Supplementary Information) which possess a wavelength 
of Xp/2 = \/2ao; with being the Fermi wavelength 
of the hybridized Fermi surface along the diagonal. The 
envelope of these oscillations decays exponentially with 
distance, r, from the Kondo hole, i.e.. Ax oc exp (— r/^). 
Such an exponential decay is expected since within the 
KBdG formalism, a Kondo hole is equivalent to an /- 
electron state with on-site energy e/ = oo. Such a state 
it localized, and its effects on Ax or As necessarily de- 
cay exponentially. With increasing magnetic interaction, 
^ increases approximately linearly (see Fig. S2 in the Sup- 
plementary Information), leading to a the growth in the 
spatial extent and amplitude of the perturbations in Ax, 
as can directly be seen from a comparison of Figs. [T^ 
and [T]:. These perturbations are therefore a direct mea- 
sure for the strength of the magnetic interaction and the 
resulting correlations in the material. 



The spatial perturbations in the electronic correla- 
tions, s(r) possess a direct spectroscopic signature in the 
local density of states of the conduction band, Nc{r,uj), 
which can be probed via scanning tunneling spectroscopy 
(STSj^SMl. The LDOS describes the quantum mechani- 
cal probability to add (for > 0) or remove (for a; < 0) 
a conduction electron with energy huj from the system 
at site r. Nc{r,uj) for an unperturbed Kondo lattice, 
shown in Fig. [2^, exhibits a gap which is the primary 
spectroscopic signature of the hybridization between the 
conduction band and /-electron state of the magnetic 
atoms, and the concomitant screening of the magnetic 
moments. This gap was recently observed for the first 
time in STS experiments by Schmidt et al.'^^. The peak 
in the LDOS at the low-energy site of the gap arises from 
the Van Hove singularity of the hybridized Fermi sea (in- 
dicated by the green arrow in Fig. [i}d). A comparison 
with the LDOS at the site of the Kondo hole reveals a 
significant shift of spectral weight from negative to posi- 
tive energies. To understand this redistribution, we note 
that the screening of a single Kondo atom by conduction 
electrons leads to an increase in the local electron density. 
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Figure 2: Spectroscopic Fingerprints a. The local con- 
duction band density of state Nc{r,uj) for the unperturbed 
Kondo lattice, and at the site of the Kondo hole, R — (0, 0), 
and its next-nearest-neighbor site R' = (1, 1) where As(R') > 
0. b. The non-local /-electron spin susceptibility, r/(r,r',aj) 
between nearest neighbor sites in the unperturbed Kondo lat- 
tice, and in the Kondo hole system between sites r = (1, 1) 
and r' = (1,0) with Ax > and between sites r = (1, 1) and 
r' = (1,2) with Ax < 0. 



71c (r), which is related to the LDOS via 

POO 

nc{r) — / dui np{u)) Nc{v,ijj) (1) 



with np being the Fermi distribution function. An in- 
crease in nc(r) therefore necessarily implies a redistribu- 
tion of spectral weight in Nc{r,oj) from positive to neg- 
ative energies. The removal of a Kondo atom from the 
lattice simply leads to the opposite effect with a decrease 
in nc(r) at the site of the hole, and the corresponding 
shift in the LDOS shown in Fig. [2^. This suggests that 



Nc{r,uj) directly reflects the perturbations in s(r), since 
for a site with As(r) > one has Anc(r) > 0, and a con- 
comitant redistribution of spectral weight in Nc{r,uj) to 
negative frequencies (and vice versa). Indeed, a plot of 
Nc{r,Lu) at the next-nearest neighbor site of the Kondo 
hole, R', with As(R') > 0, exhibits an increase in the 
spectral weight at negative frequencies, and hence con- 
firms this conclusion (see Fig. [2^) . 

The spatial perturbations in the magnetic correlations, 
x(r,r'), possess a spectroscopic signature in the non- 
local /-electron spin susceptibility, r/(r,r',w), between 
nearest neighbor sites r and r'. Tf can in general be 
measured via atomic force microscopy'^*'', and its Fourier 
transform describes how magnetic correlations decay in 
real time. A plot of Imr/(r,r',a;) in Fig. illustrates 
the important relation with Ax. In comparison to the 
unperturbed Kondo lattice, one finds that for two sites 
r,r' with Ax > 0, llmF/j is enhanced, while Ax < 
leads to its suppression. 

Replacing_ a Kondo atom by a non-magnetic 



impurity 



' gives rise to a local scattering potential for 



the conduction electrons (see Sec. S2 of the Supplemen- 
tal Information). This, in turn, induces perturbations in 
Aric, As, and Ax (see Figs. [3^-c) that possess several 
distinct differences to those caused by a Kondo hole. In 
particular, if the non-magnetic scattering potential is at- 
tractive, Uo < 0, the spatial pattern of Aric is inverted, 
i.e., a site with Arte > for the Kondo hole case, now 
has Aric < (cf. Figs. [T]i and|3^). Since the same in- 
version occurs for the spatial patterns of As and Ax, we 
conclude that the electronic and magnetic correlations 
are strongly affected by the spatial redistribution of Uc- 
Moreover, an impurity bound state can be formed around 
the non-magnetic impurity, whose spectroscopic signa- 
ture is a sharp peak in Nc{r,uj) inside the hybridization 
gap, as shown in Fig.jSji. The emergence of this impurity 
state is tied to an attractive scattering potential, Uq < 0, 
whose magnitude exceeds a critical value, Uc, which is 
determined by the particle-hole asymmetry of the con- 
duction band. A non-zero Uc could be responsible for 
the observation that different non-magnetic impurities 
can lead to disparate properties of the system^°. With 
increasing \Uo\, the bound state first emerges at the high 
energy side of the hybridization gap and then moves to 
lower energies, as shown in Fig. |3|i. This bound state 
is predominantly formed by /-electron states, as can be 
deduced from a spatial contour plot of Nc{r,ui) at the 
bound state energy, fi^, shown in Fig. [3^. The bound 
state is spatially isotropic, and decays exponentially with 
distance from the non-magnetic impurity with a decay 
length, w 0.65ao. If this bound state were formed 
by (light) conduction electron states, a simple estimate 
shows that it would possess a decay length in excess of 
60ao, in disagreement with the results of Fig. |3^. This 
result is particularly interesting since the non-magnetic 
impurity leads to scattering in the conduction band only. 
In contrast, the spatial oscillations of iVc(r,a;) for fre- 
quencies outside the hybridization gap are delocalized. 
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Figure 3: Local density of states around a defect Spatial 
contour plot of a. As, b. Ax, and c. Aric for a non-magnetic 
impurity with Uo = — l.OiJo and J/J = 0.001. d. Comparison 
of Nc{r, uu) at the site of a Kondo hole and of a non-magnetic 
impurity. An impurity bound state inside the hybridization 
gap emerges only for an attractive non-magnetic impurity if 
|C^o| > t^c = 0.62Eo. To facilitate the comparison, the curves 
of the bound state have been multiplied by 0.03. Spatial 
contour plot of A'^c(r,i^) for Uo ~ — 1.0i?o at e. the energy of 
the bound state, u = ^1b, and f. outside the hybridization 
gap at cj = 0.02i5o. 



as exemplified in Fig. |3f, and hence are predominantly 
arising from the (hght) conduction electron band. 

The strongly correlated nature of the Kondo lattice 
also manifests itself in highly non-linear quantum 
interference of the spatial perturbations in x ^-nd s 
emanating from adjacent Kondo holes. In a periodic 




Figure 4: Novel Ground State driven by Quantum In- 
terference Periodic array of Kondo holes: a. Free energy for 
an unperturbed Kondo lattice (KL) and a Kondo hole (KH) 
array as a function of 1/ J. For the electronic band structure 
considered here, one has Ic/J = 0.0011(1), indicated by an 
arrow. Contour plots of b. As, c. Ax, and (d. Aric for 
I/J = 0.0013 > Ic/J. 



array of Kondo holes, this non-linearity drives the 
system to a novel inhomogeneous ground state via a first 
order phase transition once / exceeds a critical value, I^, 
as demonstrated by a plot of the free energy in Fig. |4^. 
Note that Ic increases with increasing distance between 
the Kondo holes. This novel inhomogeneous ground 
state is characterized by spatial patterns of Ax, As, and 
Auc [see Figs. |4}3-d], that are strikingly different from 
the ones for I < Ic shown in Fig. [T] The similarity of 
these three spatial patterns, as well as the significant 
amplitude of the perturbations in the conduction elec- 
tron density, ric, suggest that while the phase transition 
is driven by quantum interference of the correlations' 
spatial perturbations, the resulting real space patterns 
are determined by the redistributed conduction electron 
density. Changing the symmetry of periodic array, for 
example, by adding a second hole to the unit cell (see 
Fig. S3 of the Supplemental Information), leads to a 
significant change in the spatial pattern for I > Ic, 
reflecting the strong non-linearity of the system. 
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Supplemental Online Information for 

Defects in Heavy- Fermion Materials: 
Unveiling Strong Correlations in Real Space 

by Jeremy Figgins and Dirk K. Morr 

Section SI: Supplemental Figures 





Figure 1: Defect in a Kondo lattice Blue dots represent the sites of the metallic conduction lattice, while red arrows show 
the locations of the magnetic moments, i.e., the magnetic Kondo atoms, in the Kondo lattice. The green ring denotes the 
site of a missing Kondo atom, i.e., a Kondo hole, or of a non-magnetic impurity, where the latter leads to scattering of the 
conduction electrons only. The defect induces spatial perturbations in the electronic and magnetic correlations of the system, 
which are determined by the real space oscillations of I^PcP with wavelength X'p/2 and with wavelength )^p/2. Here, ^'c 

and ^P/ are the wave functions of the conduction and /-electron states, respectively, with Xp and being their associated 
Fermi wavelengths. 
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Figure 2: Spacial oscillations of Ax a) Spatial oscillations of Ax(i',r') (shown as a function of Ar = |r + r'|/2 measured 
from the Kondo hole) along the lattice diagonal. The oscillations exhibit a wavelength Xp/2 and decay exponentially as e~^^^ 
(the exponential fit is shown as a red line), b) Dependence of ^ on the magnetic interaction, /. Errors in ^ are smaller than 
the data points. ^ increases approximately linearly with 1/ J (dashed line is a guide to the eye). 
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Figure 3: Quantum Interference in a periodic array of Kondo holes with two Kondo holes per unit cell Contour 
plot of (a) As and (b) Ax for I / J = 0.0013 > Jc/J. The spatial interference pattern for J > Jc is distinctively different from 
that of a periodic array with a single Kondo hole per unit cell (shown in Fig. 4), reflecting the strong non-linearity of the 
system. 
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Section S2: The Kondo-Bogoliubov-de Gennes (KBdG) Formalism 

Starting point for the study of defects in Kondo lattice systems is the Kondo Heisenberg Hamiltonian 

^ = - ^ ir,r'4^,Cr^, + J^'sf • + ^'/r,r'Sf • , (SI) 



r.r ,0" 



where try is the fermionic hopping element between sites r and r', and cl^,Cr,a creates (annihilates) a conduction 
electron with spin a at site r. J > is the Kondo coupling, and and are the S = 1/2 spin operators of 
the magnetic (Kondo) atom and the conduction electron at site r, respectively. I^y > is the antiferromagnetic 
interaction between Kondo atoms at sites r and r', which we restrict to nearest-neighbor sites below. We take the 
lattices of the conduction band and that of the magnetic (Kondo) atoms to possess identical structures and lattice 
constants. Thus, the site of a Kondo atom and the (single) conduction site that it couples to are denoted by the same 
r. The unprimed sum runs over all sites of the conduction lattice, while the primed sum runs over all positions of the 
Kondo atoms. These two sums are identical for an unperturbed Kondo lattice without any defects. 

A systematic large- iV expansiorPISl of the Hamiltonian can be achieved by generalizing the spin operators to SU{N) 
and representing them using Abrikosov pseudofermions 

Sf = ^ fla(^a,pfrj = ^ 4 „Cr„_^Cr,/j , (S2) 

a,f3 a, 13 

where a,/3 = 1,...,A'' and (Ta^p are the generators of SU{N). Here, /^^ (/r,a) creates a pseudofermion in the /- 
electron state of the magnetic Kondo atom. In addition, one needs to satisfy the constraint that each magnetic atom 
(i.e., each /-electron site) is occupied by N/2 particles, i.e., 

n/(r)-^/t„/r,. = y . (S3) 



Inserting the representations of Eq.(S2| into the Hamiltonian,Eq.(Sl ), yields quartic fermionic interaction terms. 
Since our objective is the study of Kondo lattices containing defects which break the translational invariance, it is 
necessary to introduce a mean-field decoupling of the resulting interaction terms in real space. The me an-fi eld theory 
of the homogeneous Kondo lattice (or the single Kondo impurity) becomes exact in the limit N — > ocPJ^I. However, 
it was argued that even fo r N = 2, the qualitative and to a large extent quantitative features of the mean-field 
solutions remain unchanged^, which in certain limits has been validated by the comparison with exact numerical 
renormalization group studie^^^. Following this argument, we set = 2 below, and introduce the following local and 
non-local expectation values 

= { E(/r,aCr.a) ; X(r, r') = ^(/r.a/r'.o) • (S4) 

a Q 

Here, s(r) describes the local hybridization between the conduction electron states and the magnetic /-electron 
states, whose magnitude is a measure of the Kondo screening. Thus s(r) = represents an uns creened magnetic 
moment at site r. x('';i"') is a measure of the magnetic correlations between Kondo atom^^. To enforce th e 



constraint (fi/(r)) = 1 within our formalism, we add the term X^r a ^/('')/r a/r,a to the Hamiltonian in Eq.(Sl|, 
where e/(r) represents the on-site energy of the /-electrons. The resulting Hamiltonian is quadratic, containing the 
set of parameters {s(r), x(r, r'), £:/(r)}, and therefore can be diagonalized in real space (we assume periodic boundary 



conditions). After each diagonalization, the mean-fields s(r) and x(r,r') are computed self-consistently via Eq.(S4|, 
and e/(r) is chosen such that (n/(r)) = 1 at each site. This procedure is repeated until a self-consistent solution for 
{s(r), x(r, r'), ej(r)} is obtained. We note that for an unperturbed Kondo lattice, this formalism is identical to the 
saddle-point approximation of the path-integral approactl^. 

In order to describe a Kondo hole at site R, we remove the corresponding spin operator, from the Hamiltonian 



in Eq.(Sl W- When a Kondo atom at site R is replaced by a non-magnetic impurity, the spin operator is removed, 
and the term Uq a'-R-," added to the Hamiltonian, where Uq is the non-magnetic scattering strength. Here, 

we assume that the non-magnetic atom leads to scattering of the conduction electrons only, since the magnetic 
/-electron states are not expected to hybridize with the electronic states of the non-magnetic impurity, in agreement 
with Ref. ITTl but in contrast to Ref. [12] 
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For the results shown in the article, we took the conduction and Kondo lattices to be finite two-dimensional square 
lattices containing each Af = M x M sites. Our formalism can equally well be applied to three-dimensional lattices, 
however, we do not expect that the qualitative nature of the results shown in this article will be changed. For the 
conduction band, we employ a nearest- neighbor hopping, t = 0.5i?o, and a chemical potential, fi = — 1.809i?o- Here, 
Eq is an overall energy scale related to the conduction bandwidth W = 4:Eq. These parameters yield an isotropic 
conduction Fermi surface (see Fig. lb) with Fermi- wave-length = lOao) where Uq is the lattice constant, and a 
band filling (per spin degree) of tt/IOO. 

We choose the size of the lattice sufficiently largo such that the set of solutions for {s(r), x(r, r')j ^/(i")} is 
independent of A/". The results shown in the main text were obtained for systems with M = 41, and we did not 
find any significant quantitative changes for lattice sizes up to M = 71. For the unperturbed Kondo lattices with 
J = l.OEo, the set of parameters for several values of 7/ J are shown in Table 1. 



I/J 


£/ 


s 


X 


0.0005 


0.001247 


0.048470 


0.00008122 


0.001 


0.001232 


0.048470 


0.00016579 


0.0011 


0.001229 


0.048470 


0.00018281 


0.0013 


0.001225 


0.048470 


0.00021688 



Table I: Values for s,x, and e/ (in units of Eo) for an unperturbed Kondo lattice as a function of I/J. Note that for these 
small values of I/J, s remains practically unchanged. A significant reduction of s is only seen for larger values of I/J. 



Section S3: The Local Density of States of the Conduction Band 

In general, the local density of state for the conduction band, Nc{r, oj), can be directly obtained from the eigenvectors 
and energies of the diagonalization procedure described in the previous section. However, the finite size of the systems 
leads to (artifactual) strong oscillations in the frequency dependence of the LDOS. Wc developed an alternative 
approach to circumvent this problem, and verified that it significantly reduces the artifactual frequency oscillations, 
and that it yields identical results to the first approach when considering the spatial form of Nc{r,ij) at fixed co. 
In this approach, wc compute the local conduction electron density of states at a site r for a given set of solutions, 
{s(r),x(r,r').£/W} via 

N,{r,uj) = -hm[G{v,v,uj)] , (S5) 

where G(r, r, ui) is the full retarded conduction band Green's function given by (in what follows, all Greens functions 
are retarded) 

G(r,r,w) = Go(r,r,w) + G[,''(Ar,w)si^(w)sG[,'^(-Ar,u;) . (S6) 

Here we defined 

(Go(ri,r,a;) \ 
: . (S7) 

Go(rp,r,w) / 

where Go is the Greens function of the decoupled (unhybridized) conduction band, rj is the position of the i'th Kondo 
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atom, and V is the number of Kendo atoms in the system. Moreover, for the hybridization matrix s one has 

/ s(ri) \ 

' s(r2) ' 

... 

s(rp_i) 

V s{y-p) J 



(S8) 



and F is the full /-electron Green's function with the {ij) element of its inverse given by 



uj - ef{vi) - s{vi)G(i{Y.i,Vj,uj)s{v.j) 
X(ri,rj) - s(r,)Go(r,,rj,a;)s(rj) 



\ii= j 
if i ^ j 



(S9) 



In the presence of a non-magnetic impurity scatterer at site R with potential Uq, it is easy to show that the full 
conduction electron Green's function can be obtained from 



GjvAf = G(r,r,a;) + 



G(r,R,a;)C/oG(R,r,a;) 
1 - C/oG(R,R,w) 



(SIO) 



with G(r, r'jCj) obtained from Eq.(S6). 

Section S4: The /-electron Spin Susceptibility 



The /-electron spin susceptibility, F j , in Matsubara time is given by 

r/(r,r',r) = (r.Sf(r).Sr^(0)) . (Sll) 
After Fourier transformation into Matsubara frequency space, and analytic continuation, one obtains [ks = 1) 



(S12) 



where F is the /-electron Greens functions from Eq.(S9). Using next the Lehmann representation of F^ one obtains 
for the imaginary part of the retarded /-electron spin susceptibility 

3 

ImF/(r,r',fi) = — / ImF(r, r', u) ImF(r', r, u - f]) . 
271" Jo 



(S13) 
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